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We utilize Polyakov loop correlations to study (3+l)D compact £7(1) 
flux tubes and the static electron-positron potential in lattice gauge theory. 
By using field operators it is possible in U(l) lattice gauge theory to probe 
directly the electric and magnetic fields. In order to improve the signal-to- 
noise ratio in the confinement phase, we apply the Liischer- Weiss multilevel 
algorithm. Our code is written in CUDA, and we run it in NVIDIA FERMI 
generation GPU's, in order to achieve the necessary performance for our 
computations. 



PACS numbers: 11.15.Ha; 12.38.Gc 



1. Introduction 

There are several models trying to describe confinement features that can 
achieve good successes, mainly dual superconducting models and effective 
string models. Though, the exact details of confinement mechanisms are 
still an open question in physics. 

It is important to compare the predictions of these models with more fun- 
damental calculations in order to check the domain of validity of the models 
and to understand more deeply the processes underlying confinement. 

Studies have been conducted to pursue this objective with different gauge 
groups, mainly studying the potential of two static quarks. More recently, 
some studies have also focused the predictions for the tube flux itself. Usu- 
ally this can only be done through the measurement of the squared field 
average, but in U(l), being abelian, it is possible to study the electric and 
magnetic fields directly also. 
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1.1. Effective string model 

According to the effective string model the quarks are bounded by a 
string that confines them in a asymptotically linear potential. 

This model achieved several successes, being the most well-known the 
existence of the universal 1/r term in potential (Liischer term [1]) which 
was confirmed in many studies regarding different gauge groups. 

More precisely effective string model predicts the following potential at 
leading order [1] 

V(r) = fi + ar-^=^ + 0(±) (1) 

where fi is an arbitrary constant, a is the string tension and d is the dimen- 
sion of space. 

It is now possible to check other predictions of this model including the 
ones related to the shape of the flux tube itself. 

This have been done in 2+1 dimensions with several groups, since in 2+1 
dimensions the theoretical predictions are well established, and the results 
have been found in agreement with this model. In 2+1 dimensions there 
are two results, for the limits of zero temperature and high temperature 
(although still in confining phase). These results consist of a logarithmic 
increase in the flux tube width at zero temperature and a linear increase at 
finite temperature (cf. [2]), where the flux tube is evaluated according to 
the following definition 

where 0(x±) is the profile of the field in the mediator plane of the charges. 



2. Methods 

We consider an euclidean 3+1 dimensional U(l) lattice gauge theory 
with periodic boundary conditions. The action used is Wilson action for 
U{\) 

x n<u 

where (3 = ^, U fiu (x) = U^(x) U v (x + fi) U*(x + u) U*(x) is the plaquette 
along the directions fi and v and U^x) is the gauge group element in position 
x and direction fi. 

Wilson action is local (in the sense it only depends on plaquettes around 
it) which is a condition for using the multilevel method [3] we implement. 
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2.1. Configurations 

Our lattice is an isotropic lattice of size 24 4 and we use (3 = 1. This choice 
guarantees that we are working in confining region, at zero temperature, 
although with a (3 big enough to have relatively stable results. 

We use a combination of standard Metropolis and Overrelaxation algo- 
rithms to generate the configurations. As the time corresponding to the 
iterates between configurations is negligible compared to the time of needed 
for multilevel method we undertake more than enough iterates (100 iterates) 
in between any two used configurations to ensure those configurations are 
completely uncorrelated. 

2.2. Multilevel algorithm 

Multilevel algorithm was introduced by Liischer and Weisz [3] as a way to 
enhance the signal-to-noise ratio in large loop sizes. It consists in splitting 
the lattice in several layers where spatial links are fixed and performing 
partial updates of the rest of the lattice. 

This technique is very useful to study U(l) as the signal-to-noise ratio 
decreases extremely fast for increasing loop sizes in the confining region of 
this theory. Therefore we implement a two level multilevel algorithm [3, 4], 
allowing us to calculate correlations at a distance up to 6 lattice spacings. 

To calculate the potential we also replace the temporal direction links 
with the ones obtained from the analytical version of the multihit method 
[4, 5], which exhibit a smaller variance. 

With (P*(0)P(r)) we can calculate the static potential 

y(r) = -i-ln[(P*(0)P(r))]. (4) 

For calculating the electric and magnetic field we need to introduce a 
new operator. Electromagnetic tensor components can be probed with the 
following 

a 2 F^u = ^fplm + 0(a 3 ) (5) 

and the squared components 

a A F* v = (3 [1 - Tie P^\ + 0(a 6 ). (6) 

As we want the field in the presence of polyakov loop sources we should 
introduce them in the average 



(7) 
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where 0(x) stands for any operator we want to measure. 

To compute (P*(0)P(r) 0(x)) in the multilevel scheme we need to define 
a new operator that accounts for 0(x) 

O(x ,x,t,r) = U£(x ,t)U (x + r : t)O{x,t) (8) 

and 

TO( 2 '(i , x, t) = [T(x , t, r)O(x , x,t + l, r)\ + [O(x , x, t, r)T(x ,t + 1, r)]. 
So we have for the average field (cf. figure 1) 

(P*(0)P(r)O(x)) = i^^{ 

XQ t 

([TO^(x ,x,0)T {2 \xo,2)}...[T^(x ,N t -4)T^(x ,N t - 2)] + ... 
...+[T( 2 )(x , 0)T( 2 )(x , 2)}...[T^(x ,N t - 4)T<Q)( 2 )(xo, x, N t - 2)])} 

where V is the total number of points in the lattice. 




Figure 1: Multilevel scheme for the field. 



3. Results 

3.1. Potential 

We calculate the potential from 100 multilevel configurations, each one 
with 100 level 1 and 1000 level 2 multilevel iterates, with multihit method 
for further error reduction. We fit the results (figure 2) and extract a value 
for string tension of a = 0.16719±0. 00030, if we force the Liischer term to be 
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constant (x 2 /dof = 0.140). With Liischer term as a fit parameter we obtain 
a = 0.1666 ± 0.0022 and Liischer term 0.2743 ± 0.0605 (x 2 /dof = 0.147), 
in a good agreement with the expected value of n/12 = 0.2618. We do not 
include the first two points in the fit. 
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Figure 2: Potential 



3.2. Flux tube profile 

We calculate the E x tube profile in the middle plan between the charges. 
We use 100 multilevel configurations. The function fitted is the ansatz sug- 
gested in [2]: 



(P*P P f 



fMU/ 



A exp(— x\/s) 



1 + B exp(-x 2 ± /s) 



(9) 



(P*P) •~±i-> 1 + Dexp(-x 2 ± /s)' 

As expected we can notice that the flux tube (figure 3) gets broader at bigger 
charge separations. We quantify this result through the calculation of the 
flux tube width. 
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Figure 3: Flux tube profile 



3.3. Flux tube width 

We integrate the flux tube ansatz fits to calculate the width (cf. table 
1). The errors are estimated using a jackknife algorithm. As expected the 
flux tube width grows, although it is not possible to tell its exact form. 
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Table 1: Electric field width 
4. Conclusions 

Our results for potential are in a good agreement with the predictions of 
the effective string model. 

For flux tube profile we find an increasing width, but we are not able to 
tell the exact form of the broadening. Further research is needed, to increase 
the precision of the result with more configurations and the distance up to 
which we can compute the flux tube. 
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